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ABSTRACT 

We report conclusive verification of the detection of associated H I 21 cm absorption in 
the early-type host galaxy of the compact radio source PMN J2054^-242. We estimate an 
effective spectral line velocity width of 418 ± 20 km s _1 and observed peak optical depth 
of 2.5 ± 0.2 per cent, making this one of the broadest and weakest 21 cm absorption lines 
yet detected. For T sp j n // > 100 K the atomic neutral hydrogen column density is A/hi > 
2 x 10 21 cm~ 2 . The observed spectral line profile is redshifted by 187 ± 46 km s , with 
respect to the optical spectroscopic measurement, perhaps indicating that the H I gas is 
infalling towards the central active galactic nucleus. Our initial tentative detection would 
likely have been dismissed by visual inspection, and hence its verification here is an 
excellent test of our spectral line detection technique, currently under development in 
anticipation of future next-generation 21 cm absorption-line surveys. 

Key words: methods: data analysis - galaxies: active - galaxies: ISM - radio lines: 
galaxies 



1 INTRODUCTION 

The cold atomic neutral hydrogen content of galaxies (H I) plays 
a vital role in the evolution of galaxies, through formation of 
stars and efficient accretion of matter on to active galactic nu- 
clei (AGNs). By using the 21cm hyperfine transition of Hi, 
observed as absorption against a background radio source of 
known brightness, we can measure the distribution and kine- 
matics of line-of-sight cold gas (T ~ 100 K). Previous ob- 
servations of the 21 cm line in associated absorption often ex- 
hibit spectral line profiles with both a deep, narrow compo 
nent (~10kms~ ) and broad wings from a shallow compo 



nent (- 100km s _1 ) (e.g. |Mirabel|l989[|Morganti et al.|2001 



2005a| |Salter et al.|20T0l | Allison et al.|2012a) . From observa- 
tions at high spatial resolution there, is evidence that in some 
cases these broad, shallow components might be associated with 
a circumnuclear torus (e.g. Struve & Conway 2010, Morganti 
|et al.|20iT| >. Additionally, |Morganti et ai*j |2005b) have success- 
fully detected extremely broad absorption within gas-rich hosts, 
each of which has an existing deep, narrow absorption-line com- 
ponent. These systems are associated with high-velocity gas out- 
flows, providing a probe of the feedback between the central ra- 
dio source and the star-forming interstellar medium of its host 
galaxy. Other examples of extreme broadening include the de- 
tection of H I absorption towards the candidate binary black-hole 
4C 37. 1 1 ( |Maness et al.|2004||Morganti et al.|2009) ), with a full 
width at zero-intensity of 1600 km s~ and dominated by two 
narrow, deep peaks. However, there are currently very few de- 
tections that are dominated by broad absorption, with no (or very 
little) accompanying narrow component (see e.g. Jaffe f 1990| >. 
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In previous work, Allison et al. ( 2012a) reported the tenta- 
tive detection of broad (Awfwhm ~ 500kms~ ) Hi 21cm 
absorption towards the compact flat-spectrum radio source, 
PMN 12054^1242 ^Griffith & Wright|1993[|Healey et al.|2007> , 
using a spectral line detection method based on a Monte Carlo 
Bayesian technique. From observations using the Australia Tele- 
scope Compact Array Broadband Backend (CABB; Wilson et al. 
|2011) , and by considering the ratio of the marginal likelihoods, 
the data were found to warrant the inclusion of the spectral line 
hypothesis over the continuum-only hypothesis. Here we report 
further observations to conclusively verify the presence of this 
broad and shallow absorption line. This is one of the broadest 
and weakest 21 cm absorption lines yet detected in isolation, 
with no previous indication of absorption from an existing deep, 
narrow component. Conclusive verification of this detection pro- 
vides an excellent test of our spectral line finding method, a 
technique that can be used for detecting and parametrizing such 
spectral lines in future large-scale absorption surveys. 

Throughout this Paper we adopt a flat ACDM cosmology 
with Q M = 0.27, S1 A = 0.73 and Ho = 71 kms" 1 Mpc~\ Quan- 
tities estimated from the data are given by their mean and 1 a 
uncertainty unless otherwise stated. 



2 OBSERVATIONS AND DATA REDUCTION 

In addition to the observations in 2011 April reported by |A1- 
|lison et al.| ( |2012a| >, a further 12 h observation of PMNJ2054- 
4242, at RA(J2000) = 20 h 54 m 0r:79 and Dec.(J2000) = 
-42°42'38'.'7, was carried out in 2012 lune 17-18. The CABB 
system provides a zoom band of 64 MHz across 2048 chan- 
nels, which at the centre frequency of 1.342 GHz has a ve- 
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Figure 1. The synthesized-beam weighted spectrum, combining ATCA 21 cm data from observations in 2011 April and 2012 June. The data were 
binned to a resolution of 2 1 km s _ 1 , which is approximately three times that of the individual spectra. The optical spectroscopic redshift is indicated by 
the vertical arrow (z opt = 0.042893 ± 0.000150; [Tones et al.|2009) . 



locity resolution ~ 7kms _1 . The six-element ATCA was ar- 
ranged in the 6D East-West configuration, with baseline dis- 
tances in the range 0.08 - 5.88 km. With this configuration we 
obtained an angular scale sensitivity range of approximately 8 
- 600arcsec, and a primary beam full width at half-maximum 
(FWHM) of 35 arcmin. We observed the target source in 20 min 
scans, which were interleaved with 1.5 min observations of the 
nearby secondary calibrator, PKSB2211-388 (£i.384GHz = 
1.86 ± 0.01 Jy),nfor gain and band-pass calibration. For the cal- 
ibration of the absolute flux scale, we observed PKS B1934-638 
(Si. 384 GHz = 14.94 ± 0.01 Jy) at regular intervals separated 
by 2h. For the 2012 June epoch, we achieved a total integration 
time of 8.6 h on the target source. 

We implemented the same data reduction procedure as re- 
ported by Allis on et al.| j2012al, again subtracting the baseline 
ripple signal from other continuum sources within the field of 
view. The data were flagged, calibrated and i maged using tasks 
from the MlRIAE^Jpackage i Sault et al. 1995 I, and implemented 
using a purpose-built CABB H I data reduction pipeline. We up- 
dated our original procedure to include some minor improve- 
ments to the calibration and baseline ripple subtraction, and have 
applied these here to the data obtained both in 2011 April and 
2012 June. Assuming that the flux density of PMN J2054-4242 
has not varied significantly between the two epochs (which was 
not evident from our observations), we constructed a single, nat- 
ural weighted data cube from the combined visibility data. Since 
the ATCA does not observe at constant radial velocity (known 
as Doppler tracking), the observed frequency channels for each 
epoch correspond to slightly different velocities. To obtain ap- 
proximately uniform variance in flux across the velocity spec- 
trum (while still retaining spectral information), we binned the 
combined data to a resolution of 21 kms" 1 , which is approxi- 
mately three times that of the individual spectra. Fig.[T] shows 
the resulting synthesized-beam weighted spectrum centred on 
PMN J2054-4242. Based on the median absolute deviation from 
the median (e.g. [Whiting 2012}, the estimated per spectral- 
channel uncertainty is 1.05 mjy, which is consistent with the 
uncertainties in the spectra for each epoch and the subsequent 
velocity binning (see Table[T}. The 21cm Hi absorption can 
clearly be seen upon visual inspection, slightly redshifted with 
respect to the systemic velocity, given by optical spectroscopy 
from the 6dF Galaxy Survey (6dFGS; Jones et al. 2009|. This 
feature is seen throughout observations in both epochs, and in 
both polarizations, strongly indicating that it is not generated by 
a transient artefact such as radio frequency interference. We also 
do not see this feature in spectra that are randomly selected at 
positions away from the source. At a third epoch, during another 
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observing program in 2012 April, we obtained further low spec- 
tral resolution (230 km s -1 ) CABB data. While not useful for 
parametrizing the absorption, the 21 cm spectral line profile is 
wide enough to be seen across three channels in the data and so 
provided additional verification of our detection. 



3 ANALYSIS 

Here we provide a brief summary of the method used for 21 cm 
spectral line detection and parametrization. For a more detailed 
description please refer to work by Alliso n~et al.| ( |2012a|b| , 
and references therein. In order to verify whether 21 cm H I ab- 
sorption has been detected with significance, we estimated the 
conditional probability of a spectral line model (parametrized 
by a sum of spectral line components and a first-order polyno- 
mial) and compared it to that of a continuum-only model (also 
parametrized by a first-order polynomial). We also included a 
nuisance parameter, which represents systematic error due to im- 
perfect data reduction. The model spectrum is multiplied by this 
parameter, with a non-uniform prior probability given by a nor- 
mal distribution of /i = 1.0 and a — 0.1, which represents a 
10 per cent systematic error. Note that we have revised our esti- 
mate of the systematic error in the ATCA data with respect to 
that given by |Allison etaT] l |2012a| l. 

We have assumed no prior information about the probabil- 
ity of detection of a spectral-line, since it is difficult to draw 
strong statistical conclusions from the existing small sample of 
detected associated H I absorpt ion systems (e.g.|Curran & Whit-| 



|ing|2010||Curran et al.|20TTa||Allison et al.|2012a| >. Therefore 
the ratio of probabilites for each model hypothesis, given the 
data, is equal to the ratio of probabilities for the data given the 
model hypotheses (also known as the ratio of marginal likeli- 
hoods or Bayesian evidence). We define the statistic R by 



R = In 



Em 



(1) 



where Em and i5 CO nt are the Bayesian evidence for the spectral 
line and continuum-only models respectively. If the value of R is 
greater than zero, then this gives the level of significance for the 
spectral line hypothesis. However, if the value of R is less than 
zero, then the data do not warrant the inclusion of a spectral line 
component in the model and so we reject this hypothesis. 

In order to calculate the Bayesian evidence, we used an 
efficient Monte Carlo algorithm provided by Feroz & Hobson 
(2008) and |Feroz et al.| ( |2009| >, called MultiNest, which im- 
plements the nested sampling algorithm by Skilling (2004). We 
optimized the significance statistic R by increasing the number 
of spectral line components, thus giving the best-fitting spectral 
line model. 
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Figure 2. Results of fitting to spectra from 21 cm observations in 201 1 April and 2012 June with the ATCA. The radial velocity axis is given in the rest 
frame defined by the optical spectroscopic redshift (vertical dashed line). The data have been simultaneously fitted by a combined continuum polynomial 
and a spectral line model. The solid black line represents the data after subtraction of the best-fitting (maximum-likelihood) continuum component. The 
red line represents the best-fitting spectral line model (see text for details). The blue line represents the best-fitting residual, including a velocity-axis 
offset for clarity. 



4 RESULTS AND DISCUSSION 

4.1 Parametrization of the 21 cm absorption 

We implemented the analysis method outlined in Section[3] ini- 
tially representing the 21 cm spectral line by multiple Gaussian 
components. The data in both epochs warrant the inclusion of the 
spectral line model over the continuum-only model and, by com- 
paring the Bayesian evidence for different numbers of spectral 
line components, the single-component model is favoured. How- 
ever, while a single Gaussian does provide a good fit to the com- 
bined data (at maximum likelihood xL/d.o.f = 0.94 ± 0.06), 
upon visual inspection the spectral line profile appears more 
steep sided than predicted by this model. To test this hypothe- 
sis, we tried an empirically motivated spectral line model, which 
provides for a simultaneously broad and steep-sided trough-like 
profile, and is given by the product of two Gauss error functions 
(R. Jurek and T. Westmeier, private communication). This is a 
simplified version of the 21cm spectral line model developed 
by Westmeier, Jurek & Obreschkow (in preparation). Since this 
model contains one more parameter than a single Gaussian pro- 
file, which determines the steepness of the profile sides, it must 
provide a significantly better representation of the data in order 
to be preferred. 

Again using Bayesian inference, we found that the 2011 
April data favour the single Gaussian model, while the 2012 
June and combined data favour the trough-like model. This re- 
sult is perhaps unsurprising given the relatively low signal-to- 
noise ratio (S/N) of the 201 1 April data. In the case of the com- 
bined data, the trough-like model is found to be approximately 
35 times more likely than the single Gaussian model, with a sig- 
nificance above the continuum-only model of R — 109.5 ± 0.1. 
When we replaced the uniform prior for the redshift of the 21 cm 
spectral-line with a normal prior of 12859 ± 500 kms -1 , given 
by the 6dFGS optical redshift (Jones et al.| 2009 ) and typical 
infall/outflow gas velocities, the spectral line significance in- 
creases to R — 111.8±0.1. Fig.|2]shows the best-fitting spectral 
line model for each of the individual spectra and the velocity- 
binned combined data. 

For a sufficiently bright background source, the actual op- 
tical depth is related to the continuum flux density (S CO nt), the 
fraction of flux that is covered by the foreground H I gas (/) and 
the observed spectral-line depth (ASii nc ). This relationship is 
given by 



where t(v) is the optical depth of the absorbing gas as a function 
of the rest-frame radial velocity. If we assume that the 21cm 
absorption is optically thin (ASii nc / f Scout < 0.3), the above 
expression reduces to 



t(v) 



f Scot. 



T b S (v) 

f ' 



(3) 



where r b s («) is the observed optical depth. Since many 21 cm 
absorption lines in the literature are far from being well mod- 
elled by a single Gaussian component, for the purpose of com- 
parison, we define an effective width by 



Av eS 



/ T obs (v)dv 

7~obs,pcak 



(4) 



where r b s , pC ak is the peak observed optical depth. This defini- 
tion of the spectral line width avoids missing the contribution 
of broad, shallow components (as when using the FWHM), or 
the width parameter being strongly dependent on the spectral 
uncertainty (as when using the full width at zero-intensity). For 
the special case of a single Gaussian component, the effective 
width is just a factor of 1.06 times the FWHM. In Table[T] we 
summarize the derived spectral line parameters from the best- 
fitting models. Parameter values estimated from the 2011 April 
and 2012 June data are mostly consistent with each other and 
any significant differences are the result of the choice of model 
parametrization (the single Gaussian is intrinsically more peaked 
than the trough-like model). The small difference between the 
value of R for our re-analysis of the 2011 April data, with that 
from the previous analysis by Alli son et al.| j2012a), is the result 
of improvements to our data reduction method. 

The column density of atomic neutral hydrogen (in cm -2 ) 
as a function of the integrated optical depth (in km s _1 ) is given 
by (Wolfe & Burbidge|1975) 



Nm = 1.823 x 10 18 T spin j r(v)dv 

« 1.823 x 10 18 ^p J r ohs (v)dv, 



(5) 



where T sp i n is the mean harmonic spin temperature of the gas (in 
K). Based on this relation, we estimate the Hi column density 
towards PMN J2054-4242 to be 



~{v) = - In 1 



AS linc 

fScont / 
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(2) iV H i~ 1.91 ± 0.14 xl0 21 f^^ 



cm 



(6) 
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Table 1. A summary of the derived properties of H I absorption from model fitting to the ATCA 21 cm data, where cx c han is the estimated uncertainty 
per channel; C3p ea ic is the redshift at peak spectral-line depth; ASii ne!pea k is the peak spectral line depth; S con t,peak is the continuum flux at peak 
spectral line depth; J r i, s (v)dv is the velocity integrated observed optical depth; Av c ff is the effective velocity width (as defined by equation|4j; 
R is the natural logarithm of the ratio of probabilities for the single-component spectral line model versus the continuum-only model; x^/d.o.f is 
the reduced chi-squared statistic for the best-fitting (maximum likelihood) model parameters (under the assumption of model linearity). Note that the 
uncertainties include the systematic error in the data represented by a nuisance parameter with a normal prior. 



Epoch 


Resolution 


ff chan 


^pcak 


^^linejpeak 


•Scent, peak 


jT obs (v)dv 




R 


xL/d-o.f 




(kms" 1 ) 


(mjy) 


(kms" 1 ) 


(mJy) 


(mJy) 


(kms" 1 ) 


(kms" 1 ) 






2011 April 


7 


3.83 


13073 ± 38 


4.8 ±0.8 


160 ± 11 


15.8 ±2.8 


540 ± 120 


24.2 ±0.1 


0.98 ±0.03 


2012 June 


7 


1.91 


13048 ± 10 


3.9 ±0.4 


163 ± 10 


10.1 ± 0.7 


418 ± 20 


95.0 ±0.1 


1.02 ± 0.03 


Combined 


21 


1.05 


13046 ± 10 


4.1 ± 0.4 


163 ± 10 


10.5 ± 0.7 


418 ± 20 


109.5 ±0.1 


0.91 ± 0.06 
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^cm" 2 ). 



4.2 High velocity cold gas towards PMN J2054^1242 

By applying Bayesian inference to our 21cm ATCA data we 
have detected cold, high column density Hi gas towards the 
compact flat-spectrum radio source PMN J2054—4242. The ob- 
served total integrated optical depth for the absorption is rel- 
atively large (see Table[TJ, however, this is distributed over a 
broad spectral line width, producing a modest peak optical depth 
of T bs,pcak = 2.5 ± 0.2 percent. The 21 cm spectral line is red- 
shifted with respect to t he optical spectr oscopic redshift (z op t 
0.042893 ± 0.000150; 



Jones et al.|2009| by 187 ± 46 km s 



possibly indicating that the gas is infalling towards the central 
radio source. The strong broadening of the line perhaps sug- 
gests that the cold gas is rotating with very high velocity, or 
that the infall is accelerated along the line of sight. A high 
spatial resolution 2 1 cm study of the similarly compact radio 
source PKSB1814-637, by |Morganti et al.| ( |20li) , showed that 
the broadened absorption component in that system was likely 
located within a circumnuclear disc, while the deeper narrow 
component arises from absorption in an extended disc of gas on 
the kpc scale. The absence of a narrow component in the spec- 
trum of PMN J2054^-242 is perhaps indicative of either low 
column density (or missing) large-scale cold gas, or that any 
extended structure is orientated away from the line of sight to 
the radio source. These hypotheses are consistent with the opti- 
cal classifications of the host galaxies for each of these sources: 
PKSB 1814-637 is hosted by a gas-rich regular edge-on disc, 
while for PMN J2054-4242 the host is an early type ( |Loveday| 
1996). However, without further information, we can only spec- 
ulate here as to the physical nature of the spectral line broaden- 
ing. We are therefore pursuing follow-up multiwavelength and 
high spatial resolution observations to test these interpretations 
of the 2 lem ATCA data. 



4.3 The scarcity of broad, weak 21 cm absorption lines 

The left-hand panel of Fig.|2]shows the initial detection of broad 
21 cm absorption towards PMN J2054-4242 (2011 April). From 
visual inspection of this spectrum, one might reject the feature 
as being an artefact, due to its wide profile and low peak S/N. 
Indeed by using the automated source-finding tool DUCHAMP 
jWhiting|2012) , this original feature was only detected with sig- 
nificance once we had smoothed the spectral data (using a Hann 



3 The lowest value of T sp i n // found so far, for a damped Lyman-a 
absorber, is 60 K jCurran et al.|2007| . 



window function with a width of 5) and applied a threshold cut- 
off of three times the S/N. However, by using a Bayesian ap- 
proach to spectral line finding and fitting, A llison et al.| ( |2012"a) 
found that the feature was significant above the noise and that it 
occupied a parameter space that was different to that of the false 
detections from continuum artefacts. The redshift is close to that 
of the optical spectroscopic measurement, which strongly sug- 
gests that it arises from associated H I absorption. All of these 
indicators prompted further observations in 2012 June, with a 
longer integration time, leading to confirmation that the feature 
is indeed a 21 cm absorpion line. 

Such broad, shallow 21 cm absorption lines physically cor- 
respond to a large spread in the velocity of H I gas along the 
line of sight. Existing detections of broad absorption are often 
associated with an acompanying narrow, deeper component. For 
example, Morganti et al. (2005b) conducted a targeted survey 
of 11 bright radio-loud AGNs, with evidence of past star for- 
mation and a rich inter-stellar medium, and most with existing 
detections of deep, narrow H I absorption. The broad bandwidth 
available (±2000 kms" 1 ) and the high sensitivity achieved (rms 
~0.4 - 0.8 mJy beam" 1 ), using the Westerbork Synthesis Ra- 
dio Telescope, allowed detection of broad and shallower absorp- 
tion with Tp Ca k < 0.01. Morganti et al. detected broad, shal- 
low absorption towards 6 radio galaxies, all of which were as- 
sociated with existing deep, narrow components, and had bulk 
blueshifted motion interpreted as line-of-sight cold gas outflows. 
However, our detection of a single broad absorption line towards 
PMN J2054-4242 raises the question as to whether there exist a 
population of very shallow and broad H I absorption lines, which 
have no associated narrow, deep component, and as yet remain 
largely undetected. 

To address this question, we show in the top panel of Fig.[3] 
the peak observed optical depth, versus the effective spectral line 
width, for a representative sample of d etections of 21 cm a sso- 
ciated H I absorption in the literatur^] ( jCurran et al. |2013| also 
present a similar plot of the FWHM for a literature sample of 
redshifted absorbers). We have treated multiple spectral com- 
ponents as a single detection, so that a profile exhibiting a sin- 
gle broad component will have a larger effective width than that 
with both broad and strong narrow components. Upon inspec- 
tion of Fig.[5] it is evident that these data exhibit a weak anticor- 



4 References for the literature sample shown in F ig[3| |Shostak et "all 
(1983), |van Gorkom et al.| (T986] [T989) , |MirabeT] (1989| |1990), |Jaffe 
Uson et al.H199iy|Carilli et al.|(1992||1998) , |Morganti et al. 

2000), |Morganti et al 

2003 1, Maness 



(1990 

( 1998), |Moore et al.|(T999),|Peck et aT|(T999l 
(2001 1, Ishwara-Chandra et al. (2003 1, Vermeulen et al 
et al. (20041. Morganti et al. (2005a I, Gupta & Saikia (2006 1, Curran 
et al.| (2006), |Orienti et al. 1(2006), |Gupta et al. 1(2006), |Emonts et al. 



(2008) 7 TStruve & Conwayf (2010), (Salter et al.|p01*0) , 

ei ai'iprnV 



Emonts et al. 



(20T0) , [Curi-an et al.|(2011a|b) , |Chandoia 
J20TT 



and 



Allison et al. 



|2012a| . 



Emonts et al. 
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Figure 3. The observed 21 cm peak optical depth (top), and in flux units 
(bottom), versus the spectral line width (as defined by equation|4j, for 
detections of associated H I absorption. The empty stars represent a sam- 
ple from the literature (see text for details and references therein) and 
the filled star represents our detection towards PMN J2054^1242. The 
lines represent H I absorption for the maximum observed column density 
~ 10 22 cm" 2 , and constant T spin /f of 100K (solid), 300K (dashed) 
and 1000 K (dot-dashed). 



relation between the peak observed optical depth and spectral 
line width (the Pearson product moment correlation coefficient 
r « —0.3), consistent with the maximum column density limit 
~ 10 22 cm -2 , seen in observations of Hi gas both at z w 
jZwaan et al.|2005| > and z w 3 ( |Altay et al.|2011J . However, it 
should be noted that this figure is purely indicative of the overall 
trend, since we have not accounted for the myriad uncertainties 
in these individual estimates. 

The large majority of detections in the literature sample 
have spectral line widths smaller than our detection towards 
PMN J2054-4242 making this one of the broadest 21cm ab- 
sorption lines yet detected. The widest spectral line shown in 
Fig.[3]arises from absorption in the central gala xy NGC 1275 of 



the Perseus Cluster (~ 550km s _i ; Jaffe 1990|, which similarly 
has a profile dominated by broad absorption. However, the back- 
ground radio source in this detection (3C 84) has a flux density of 
19.7 Jy, making it a substantially more powerful probe of low op- 
tical depth gas than PMN J2054-4242. Therefore, in the bottom 
panel of Fig.[5] we have plotted the product of the peak observed 
optical depth and continuum flux density versus the spectral line 
width. Our detection occupies a region of parameter space with 
low spectral line depth and broad width, which has not yet been 
largely explored by existing searches for associated H I absorp- 
tion. Finding such a rarely detected system, in data with a rela- 
tively low S/N, demonstrates the power of our analysis method. 
The development of such sophisticated automated spectral line 
finding techniques will be vitally important for detection of these 
broad, weak absorption lines in future large-scale 2 1 cm absorp- 
tion surveys. 
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